Observation of vortex lattice melting in large untwinned YBa2Cu3C>7 _ z single crystals 
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We present a study of the vortex lattice in untwinned YBa2Cu307_i, crystals, using a combination 
of muon spin rotation and neutron small angle scattering measurements. Both methods show a very 
sharp melting temperature consistent with a first order transition. The dependence of the melting 
temperature on the angle of the field with respect to the crystallographic c-axis is studied. The 
results are compared to thermal measurements. 
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I. INTRODUCTION 



Due to their short coherence length, long penetra- 
tion depth and high anisotropy, the high-T c supercon- 
ductors (HTS), show very unusual behaviour. Apart 
from a proposed vortex-glass phase 1 , which has probably 
been observed in the highly anisotropic P^S^CaC^Os 
(BSCCO) at high fields 2 , there is also the exciting phe- 
nomenon of vortex-lattice melting in these compounds. 
Unlike in BSCCO, the vortex lattice melting transition 
in YBa2Cu307_a; (YBCO) takes place at sufficiently high 
fields to be studied using thermal measurements such as 
differential thermal analysis (DTA) 3 and specific heat 4 . 
This is mainly due to the much lower anisotropy of YBCO 
(7 = Xc/Xab — 5). Apart from this phenomenological dif- 
ference, there are however also more fundamental differ- 
ences between YBCO and BSCCO in the crystallographic 
symmetry, which is orthorhombic in the former, while it 
is tetragonal in the latter. This leads to a fundamental 
difference in the growth of the samples, in that untreated 
YBCO samples are usually heavily twinned (having do- 
mains of interchanged a and b axes). The boundaries 
between these domains (twin-planes) act as very strong 
pinning sites to the vortices 5 , which may change some of 



the properties of the vortex lattice in the sample. It has 
for instance been argued that such strong pinning may 
change the (intrinsically) first order transition to a sec- 
ond order transition 1 . This was also supported by the 
fact that the transitions observed (for instance by muon- 
spin rotation (/xSR) 2 and neutron small-angle scattering 
(SANS) 6 ) in BSCCO were generally very sharp, not in- 
consistent with a first order transition, whereas twinned 
YBCO samples gave a more gradual transition. With 
the advent of sufficiently sized untwinned and high pu- 
rity twinned samples, it became possible to study these 
effects in more detail using DTA 3 and specific heat 7 
to determine the true thermodynamic properties of the 
phase transition. In untwinned samples, DTA measure- 
ments have clearly shown the existence of a first order 
transition 3 , even though the sample wasn't pure but had 
some Au impurities due to the growth process. In ex- 
tremely pure samples, specific heat measurements have 
also shown a first order transition even though the sam- 
ple was heavily twinned 7 . However, the twinned case is 
in general much more complicated than the untwinned 
one in that the order of the transition may also depend 
on the applied field 7 . The order of the transition is thus 
probably influenced by both kinds of defects (impurities 
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and twin planes), where an excess of pinning sites of any 
kind drives the transition to second order. 

Of further interest is the dependence of the melting 
transition on the angle of the field with respect to the 
c-axis 8 . There are detailed predictions on the field de- 
pendence of the melting temperature (melting line) on 
the angle from anisotropic London theory 9 . The nature 
of the transition, especially the change in entropy at melt- 
ing has however only recently been investigated 8,10 . 

We have now performed fiSK and SANS measurements 
of the vortex lattice in a large untwinned crystal of YBCO 
with the field at different angles to the crystallographic 
c-axis. As is the case in the DTA and specific heat mea- 
surements, we probe the equilibrium state of the vortex 
lattice. From SANS, we obtain a microscopic measure- 
ment of the global (long range) arrangement of the flux 
lines. Appropriately scaled, the SANS intensity presents 
a phcnomcnological order parameter of the flux lattice 
(sec below). This allows investigations similar to those 
of specific heat or DTA, clarifying the order of the tran- 
sitions. It should however be noted that the samples 
used in SANS are very much bigger than those used for 
thermal measurements. Thus a possibly present first or- 
der transition may already be smeared by slight inhomo- 
geneities over the big sample volume. In contrast, with 
^iSR we observe the local magnetic field structure aver- 
aged over the sample. Therefore smearing effects may 
still be present, however a sharp change of the local equi- 
librium structure in itself may not be regarded as a sign 
of a first order transition. Due to the different time scales 
available to neutrons and muons one would expect to ob- 
serve a difference in the melting line determined by the 
two methods if the transition is of a more glassy nature. 
For a (first order) phase transition, one would expect 
both methods of measurement to agree. 

As will be shown below, our results differ quite strongly 
from those on a heavily twinned sample 11 ' 12 . There the 
melting transition observed with SANS is of second or- 
der, with the intensity decreasing to zero continuously at 
a temperature T m < T c . In contrast our results show 
a sharp decrease of the SANS intensity with an infinite 
slope at T m . This is consistent with a first order transi- 
tion considering the large size of the sample (see below) . 

II. EXPERIMENTAL DETAILS 

The sample consisted of a large (1.125 g) YBCO sin- 
gle crystal. It was grown untwinned by application af 
uniaxial stress during the cooling from the melt 13 . For 
the oxygenation process, the sample was treated in flow- 
ing O2 at 490°C. Magnetisation measurements using a 
vibrating sample magnetometer (VSM) indicate a T c of 
93 K. The superconducting transition is estimated from 
these measurements to be ~0.5 K. This also gives a limit 
for the sharpness of the melting transition (see below). 

The ^SR experiments were carried out on beam- 
line 7rM3 at the Paul Scherrer Institute (PSI), Villigen 
Switzerland, with a maximum field of 0.6 T and a cryo- 



stat capable of cooling to 2 K. In the experiments, a 
magnetic field was applied perpendicular to the incom- 
ing muon spin polarisation, while at a specified angle to 
the crystallographic c-direction of the sample. In this 
case, the spin of a stopped muon precesses with the fre- 
quency lu — 7 AI B/ oc , where 7^ = 2ir 135 MHz/T is the 
gyromagnetic ratio of the muon in the transverse local 
field Bi oc . The angular anisotropy of the muon decay, 
leading to a preferred emission of the decay-positron in 
the direction of the muon spin due to parity violation, 
allows the time evolution of the muon polarisation to 
be measured. In the mixed state of a type-II supercon- 
ductor, the magnetic field distribution is given by the 
structural properties of the vortex lattice (i.e. the pen- 
etration depth A, the coherence length £ and the lattice 
morphology). The vortex lattice spacing is much bigger 
than the atomic spacing of the underlying crystal struc- 
ture and is incommensurate to it. Therefore the muons 
sample the field distribution randomly, even though they 
are stopped at specific crystallographic sites. This allows 
the field probability distribution p(B) to be measured via 
the Fourier transform of the oscillating time signal. As 
the different local fields cause the muon spin to precess 
at different frequencies the time signal consists of a sum 
of oscillations at different frequencies weighted by their 
relative occurrence. In order to obtain the Fourier trans- 
form from the observed time data, we use a Maximum 
Entropy algorithm 14 . In this procedure aberrations of 
the spectra due to finite time windows or large statistical 
errors at long times are avoided, as the data are not sub- 
jected to a Fourier transform themselves. Starting from 
an initially flat field distribution p(B), the most uniform 
distribution compatible with the given time spectrum is 
obtained by maximising its constrained entropy 

L = S-\ X 2 . (1) 

Here S = ^p(B)ln(p) is the entropy of the distribution 
according to Shannon and \ 2 = ^2{Ui ~ di) 2 jo 2 is the 
usual measure giving the difference of the Fourier trans- 
form of p(B) (yi) to the data in the time domain (di), 
where a is the error in the data. 

The field distribution corresponding to a vortex lat- 
tice has very characteristic features arising from its two- 
dimensional periodicity. At extremal points in the spatial 
distribution the field probability distribution exhibits van 
Hove singularities. These singularities ideally are finite 
jumps at the minimum and maximum fields, and further- 
more a logarithmic divergence at saddle points. Due to 
the very high fields in the vortex cores, the field distri- 
bution is highly asymmetric with a long 'tail' at fields 
higher than the average. In the entangled vortex liquid 
state however, the core fields are smeared out over a large 
distance. Therefore the high field tail is truncated leading 
to a symmetric lineshape. We quantify the asymmetry 
of the lineshape by the dimensionless parameter 

(A5 3 ) 1 / 3 

a= (AB 2 ) 1 / 2 ' (2) 



A symmetric distribution corresponds to a value of a = 0, 
whereas a distribution with a high field tail has a positive 
a. 

The SANS experiments were carried out on instru- 
ments Dll, D22 and D17 at the Institut Laue-Langevin 
(ILL), Grenoble France. The neutrons typically had a 
wavelength of X n = 1 - 1.5 nm and the beam was colli- 
mated over a distance of 2.75 m in the experiments car- 
ried out at D17, whereas the collimation length was 10 
m in those at Dll and D22. The distance between the 
sample and the detector was then varied for optimal res- 
olution. The incoming neutron beam was aligned with 
the applied field to an accuracy of 0.1° by observing the 
diffraction pattern from a Nb single crystal. The sample 
was mounted in a cryomagnet capable of a field of 5 T 
and a temperature of 1.5 K. 

The neutrons are scattered from the spatial variations 
of the magnetic field in the sample, due to the magnetic 
moment of the neutron [ji n = — 1.913/xjv), with the in- 
teraction energy jj n - B. In the case of a type-II supercon- 
ductor in the mixed state, the field variation is periodic 
giving rise to Bragg diffraction off the scattering planes 
of the vortex lattice 15 ' 16 , where the angle of refraction is 
given by 

A„ = 2d hk sin(9). (3) 

Here dhk is the lattice spacing, which depends on the 
applied field as dhk — (v / 3/2$o/S) 1 / 2 for a hexagonal 
lattice ($o = h/2c is the magnetic flux quantum). Thus 
the scattering is at very small angle, as typical vortex 
separations of ^50 nm (corresponding to a field of ~1 
T) result in a Bragg angle (from Eq. 3) of 29 ~ 2°. The 
intensity of neutrons scattered into a Bragg-peak, inte- 
grated over angle as the sample is rocked through the 
Bragg condition, is then mainly given by the square of 
the form factor of a flux line, which in the London ap- 
proximation is given by 

where A is the penetration depth. The London approxi- 
mation is valid in the case of an extreme type-II super- 
conductor with A » £, a condition well satisfied by the 
HTS. Therefore, for fields much greater than B Cl , corre- 
sponding to Xq » 1, the scattered intensity depends on 
the penetration depth as / cx A~ 4 . For a classical type- 
II superconductor, the temperature dependence of the 
scattered intensity therefore is a direct measure of that 
of the penetration depth. In the HTS however, where 
thermal fluctuations may have a dominant contribution 
this is not necessarily the case. In order to discern effects 
arising from the temperature dependence of A and those 
of the thermal fluctuations finally leading to the melting 
of the lattice we divide the observed intensities by the 
temperature dependence of A~ 4 . This was obtained from 
the SANS intensity at a low field (0.2 T), where the ef- 
fects of thermal fluctuations are negligible (see below). 
This quantity then takes the role of an effective flux line 



lattice (FLL) order parameter, which we use to observe 
the melting transition and its order. 

Because the scattered intensity depends on the pen- 
etration depth as I cx A~ 4 the intensity in the HTS 
is rather low, due to their long penetration depths (A 
= 150 - 200 nm). Furthermore, substantial background 
scattering from the cryostat and extended defects in the 
sample make it necessary to subtract a diffraction pat- 
tern taken above T c to obtain a clear signal. Due to 
thermal contractions of the sample stick and cryostat, 
this background is shifted by fractions of a detector pixel 
in order to obtain a good subtraction. Thanks to the 
high neutron intensity at the instruments D22 and Dll, 
a counting time of 30 minutes for each foreground and 
background was enough to obtain a reasonable signal. 



III. RESULTS AND DISCUSSION 

A typical diffraction pattern can be seen in Fig. 1 for 
an applied field of 0.2 T. There the field is applied at an 
angle of 33° to the c-axis. This is why the observed pat- 
tern is isotropic. The distortion due to the afr-anisotropy 
(lab = 1.16 17,18 ) is compensated by the distortion due 
to the ac-anisotropy at an angle of 33°. Furthermore, 
we note that second order diffraction spots are visible 
in the figure, indicating the high quality of the flux lat- 
tice. The intensities measured from such a pattern arc 
obtained from fitting a gaussian to each individual Bragg- 
spot and taking the sum of all spots in order to reduce 
the statistical error. All of these measurements of scat- 
tered intensity rely on the fact that the width of the 
rocking curve does not change as a function of temper- 
ature. This is because the calculation relating the scat- 
tered intensity to the penetration depth is done for the 
total integrated intensity 19 , which is given by the width 
of the rocking curve multiplied by the intensity when ful- 
filling the Bragg condition. The rocking curve width was 
checked to be constant for selected temperatures over a 
wide range. We therefore only consider the scattered in- 
tensities at the Bragg-angle. 

In order to determine an effective FLL order param- 
eter we have to divide the intensity by its temperature 
dependence intrinsically given by the penetration depth. 
Therefore we have measured the scattered intensity at a 
low field (0.2 T) where the influence of thermal fluctua- 
tions should be small. The result can be seen in Fig. 2, 
where the intensity at 0.2 T is given divided by the ex- 
pectation of the 3D XY model. As the figure shows, 
the approximation of the 3D XY model, with A oc (1 
- T/T c ) -0 - 33 is valid for temperatures well below T c in 
all of the region of interest. In the determination of the 
order parameter, the intensity was thus divided by (1 - 
T/T,,) 1 - 33 for simplicity. Furthermore the insert of Fig. 2 
shows the temperature dependence over the full range. 
There we have plotted the square root of the scattered 
intensity to present the familiar A~ 2 dependence 20 . It 
can be seen that at lower temperatures the dependence 
is linear down to temperatures of ^4 K, below which it 



saturates. For a d-wave superconductor, as the HTS are 
supposed to be, a linear temperature dependence would 
be expected down to the lowest temperatures from sim- 
ple London theory. However when non-local interactions 
are taken into account 21 , a saturation at very low tem- 
peratures is predicted, not inconsistent with our data. 
We will come back to the temperature dependence of A 
below in the discussion of the melting line. However, 
we note that such a big range of applicability (down to 
^20 K below T c ) for the 3D XY model is surprising, as 
this is thought only to be valid in the region of critical 
fluctuations in close vicinity of T c . 

In Fig. 3, we present the temperature dependence of 
the FLL order parameter in the vicinity of T c for two 
different fields. At both fields, the order parameter is 
constant at low temperature and decreases gradually to- 
wards the melting temperature. When the order param- 
eter finally drops to zero, the change is very sharp oc- 
curring over a temperature range of less than 1 K. The 
gradual decrease of the order parameter can be modeled 
by including a Debye- Waller factor into the theoretical 
description of the scattered intensity. The final 'jump' 
however, is a sign of the transition being of first order. 
The width of the transition of 1 K is not much bigger 
than that observed at T c , which may be due to sam- 
ple inhomogencities. Such inhomogeneities leading to a 
distribution of T c may also lead to a distribution of T m 
which would probably be of the same order of magnitude. 
Furthermore, the distribution of fields inside the sample 
is inhomogenious 22 , as we will also see in the discussion 
of the /xSR results below, which will also lead to a fur- 
ther smearing of the melting temperatures. However, as 
this macroscopic field distribution is very narrow com- 
pared to the applied fields, this shouldn't be a big factor. 
It is therefore not unreasonable to expect the melting 
transition in such a big sample to be distributed over an 
interval of ~1 K. Note also that the order parameter (the 
scattered intensity) decreases to zero above the melting 
transition, in agreement with the expectation for an en- 
tangled vortex liquid 23,24 made up of two-dimensional 
'pancake' vortices. 

The situation for the untwinncd sample has to be com- 
pared to that in a heavily twinned sample 11 . This is 
done in Fig. 4, where we show the FLL order parame- 
ter for both a twinned and our untwinned sample in a 
field of 4 T. As can be seen in the figure, the transition 
is much sharper in the untwinned one. The order pa- 
rameter decreases more or less linearly over a range of 
5 K in the twinned sample, as would be expected from 
mean field theory. In the untwinned sample in contrast, 
the transition is at most 1 K wide. Takeing also into ac- 
count the temperature dependence at lower temperature, 
where the untwinncd sample varies much slower than the 
twinned one, the behaviour of the untwinned sample is 
consistent with a first order transition. The width of the 
observed transition is furthermore not very much bigger 
than those determined from magnetisation 25 and thermal 
measurements 8 , which are done on much smaller samples 
where inhomogeneities and possible residual twin-plancs 



must play a much smaller role. 

We now turn to the discussion of the /iSR results. Typ- 
ical field probability distributions for temperatures above 
and below the melting transition can be seen in Fig. 5 
for a field of 0.6 T. The situation is similar to that of 
BSCCO 2 ' 26 , with a lincshape highly asymmetric towards 
high fields below the melting temperature. As was dis- 
cussed above, the lineshape for an ideal lattice made up 
of straight vortex lines is expected to show such a high 
field tail. Above the melting temperature however, the 
field distribution is skewed the opposite way. This has 
been argued to arise from a macroscopic distribution of 
the local fields due to the sharp edges of the sample, with 
fields decreasing at the edges 27 ' 22 . This kind of macro- 
scopic distribution is also present below the melting tem- 
perature, however the large width of the field distribution 
due to the well formed lattice overcomes this effect, such 
that an almost ideal distribution is observed. Above the 
melting temperature in contrast, the field distribution is 
smeared by the fast thermal fluctuations of the pancake 
vortices. This leads to a strongly decreases the width of 
the microscopic field distribution ('motional narrowing') 
and hence the observed distribution follows the macro- 
scopic one. 

In Fig. 6, we present the temperature dependence of 
the /iSR lineshape asymmetry parameter a, together 
with that of the mean internal field. The mean field 
allows a direct measurement of T c via the onset of dia- 
magnetic behaviour. As can be seen in the figure even at 
the very low field of 0.3 T, we observe a sharp melting 
transition in the temperature dependence of a somewhat 
below T c . The change in lineshape displayed in Fig. 5 
manifests itself by a change in sign of a. The transition 
is still very sharp, although thermal measurements arc 
not sensitive enough to resolve a first order transition 
at such low fields. There have however recently been 
magnetisation measurements by means of torque magne- 
tometry, resolving a step in magnetisation and hence a 
first order transition at similarly low fields 28 . 

The field dependence of the melting temperatures thus 
determined (the melting line) can be seen in Fig. 7. In 
that figure we have also included the line of onset of irre- 
versible magnetisation (irreversibility line, IL) as deter- 
mined with a VSM. The microscopic melting line and the 
IL are in good agreement with each other. They follow 
a power law dependence of the form 

B m (T) = B (l-T/T c ) n . (5) 

A least squares fit to both sets of data simultaneously 
yields the values of B = 100(10) T and n = 1.30(5). This 
is in excellent agreement with many other studies 3 ' 12 ' 11 
on both twinned and untwinned samples. It is thus only 
the order of the transition, which is affected by the pres- 
ence of disorder (e.g. in the form of twin-planes) and not 
the meltinglinc itself. The shape of the melting line has 
been the ground of some debate, as a dependence Bvg ^ 
(1 - T/Tc) 1 - 33 is expected for the vortex-glass transition 
proposed by Fisher et al. 1 . This would indicate that the 
observed melting line is more likely to be a glass transi- 
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tion, which would however be in contradiction with the 
observation of a first order transition. The temperature 
dependence of a melting line is usually calculated using 
the Lindcmann criterion that melting occurs when the 
mean square displacements of flux lines make up a cer- 
tain fraction of the intervortex separation (u 2 ) = c\a\. 
The Lindcmann number cl is supposed to be a universal 
constant in all kinds of melting processes and theoretical 
estimates indicate it to be 0.1 - 0.2. Using this criterion 
and modeling the thermal fluctuations in the vortex lat- 
tice one can then derive an expression for the melting 
line 29 



m (1.57r M0 7fcBTA 2 (T))2- v ' 

Due to the mean field assumption of the temperature 
dependence of the penetration depth A~ 2 (T) oc (1 - 
T/Tc) 1 ' 2 this expression is usually taken to result in a 
value of n = 2. As we have seen however in Fig. 2, A fol- 
lows a powerlaw, such that the value of n = 1.33 is recov- 
ered over the range of temperatures where the mcltinglinc 
can be determined. Furthermore, this expression allows 
a determination of Bq, where cl is the only adjustable 
parameter. From the Intensity in the Bragg peaks, we 
have an independent determination of A = 138 nm 17 and 
the anisotropy of our sample has been determined from 
the angular dependence of the [iSR linewidth as well as 
from the distortion of the SANS pattern at high angles 
to the c-direction. Therefore our value of B = 100(10) 
T corresponds to a Lindemann number of cl = 0.173(7). 
This is in excellent agreement with that found from the 
melting line in BSCCO using ^iSR 26 . We would like to 
note at this point that the Lindemann number deter- 
mined in BSCCO relied on a different expression for the 
thermal fluctuations taking into account the electromag- 
netic coupling between pancake vortices. This indicates 
that the Lindemann criterion may indeed be universally 
applicable. 

Finally, we discuss the dependence of the melting tran- 
sition on the angle of the applied field to the c-direction. 
From the anisotropic scaling theory of Blatter et al. 9 , a 
dependence of the melting line on angle is 

B m (#) = B m (0)(cos 2 (0)+l- 2 Mn 2 m- 1/2 - (7) 

In the region of angles we investigated ($ < 50°), the 
influence of 7 on B m is very small such that it cannot be 
accurately determined this way. However using the value 
obtained from ^iSR and the distortion of the pattern 18 , 
we may calculate the expected melting lines at an angle 
from that determined parallel to c. This has been done 
in Fig. 8, where we show the melting temperatures de- 
termined by SANS at the three different angles of 0°, 33° 
and 55°. The lines correspond to the fit to the melting 
line discussed above, where the dependence on the angle 
has been taken into account as well. As can be seen in the 
figure, the data are in good agreement with the expecta- 
tions of the anisotropic scaling theory. This is similar to 
the results of DTA measurements of the melting transi- 
tion as a function of angle 8 . In addition, that work also 



provided a determination of the entropy change at the 
transition as a function of angle. Due to a cancellation 
of the angular dependencies in the melting line and the 
magnetisation jump, the entropy-change turns out to be 
independent of the angle. This is not inconsistent with 
the SANS result displayed in Fig. 9, where the tempera- 
ture dependence of the FLL order parameter is shown at 
the melting transition for different angles. The temper- 
atures are normalised to T m for direct comparison and 
the applied field was 5 T. The decrease of the order pa- 
rameter can be seen to be independent of the angle as it 
should be for a constant entropy change at the transition. 



IV. CONCLUSIONS 



In conclusion, we have presented SANS and /xSR mea- 
surements of the melting of the vortex lattice in a sample 
of untwinned YBCO. The observed behaviour is consis- 
tent with DTA measurements of the latent heat showing 
the existence of a first order transition 3 ' 8 . This is inferred 
from a phenomcnological FLL order parameter in SANS, 
which decreases to zero sharply at the melting transition. 
These findings are compared to a previous SANS study 
of melting in twinned samples 11 . The same FLL order 
parameter vanishes continuously in that case, consistent 
with a second order phase transition. 

We have furthermore studied the effects of anisotropy 
by applying the field at an angle to the c-direction. 
The melting temperatures at constant field are consistent 
with the expectations of the scaling theory of Blatter et 
al. 9 . A determination of the anisotropy 7 however is not 
possible, as the deviations of a cos (0) dependence are 
negligible in the range of angles studied. The nature of 
the transition is unchanged by having the field at an an- 
gle, consistent with the findings of Schilling et al. 8 , who 
studied the dependence of the entropy change on angle. 

The temperature dependence of the melting line is con- 
sistent with that expected from a Lindcmann criterion 29 , 
where the temperature dependence of the penetration 
depth follows that of the 3D XY model. The observed 
melting line extends down to very low fields (0.2 T), 
where /iSR still indicates a very sharp transition. The 
range of applicability of this temperature dependence is 
limited, as we have shown from the temperature depen- 
dence of the scattered intensity at low fields, where ther- 
mal fluctuations do not play an important role. For tem- 
peratures below 75 K the approximation becomes invalid. 
We would like to note that this is a rather unexpectedly 
big range of applicability of the critical exponent of the 
3D XY model. Using the independently measured pene- 
tration depth 17 and anisotropy 18 for our sample, we have 
also determined the Lindemann number from the melting 
line, which turns out to be Cl = 0.173(7). This is in very 
good agreement with the Lindcmann number previously 
obtained from a detailed /xSR study of the melting line 
in BSCCO 26 . 
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FIG. 1. A typical SANS diffraction pattern from the vor- 
tex lattice in YBCO. A field of 0.2 T was applied at 33° to 
the crystallographic c-axis and the sample was cooled through 
T c in the field to a temperature of 5 K where the image was 
taken. In order to remove a background of small angle scat- 
tering, a diffraction pattern taken above T c was subtracted. 
Note that the diffraction pattern is an undistorted hexagon 
due to the fact that the field is at 33°, undoing the in-plane 
anisotropy j a t,. Furthermore, second order diffraction spots 
are visible in the figure, indicating the high degree of order in 
the flux lattice. 
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FIG. 2. The temperature dependence at a low field of 0.2 
T gives a measure of the penetration depth. For clarity, the 
scattered intensity has been divided by (1 - T/Tc) 1 ' 33 , which 
would be the expected temperature dependence from the 3D 
XY model. This is thought to be applicable in the vicinity 
of T c . Our data however show that even at temperatures as 
low as 75 K the approximation is still valid. The insert shows 
the full temperature dependence of A~ 2 , indicating a linear 
variation at low temperatures. The saturation at very low 
temperatures may be due to non-local effects expected in a 
d-wave superconductor 21 (see the text). 
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FIG. 3. The temperature dependence of the FLL order pa- 
rameter at two different fields of 1.5 and 5 T. In both cases, the 
order parameter decreases to zero very rapidly at the melting 
transition. This is consistent with the transition being of first 
order (see the text). Furthermore the fact that the scattered 
intensity is zero above the melting temperature indicates the 
presence of an entangled vortex liquid. 
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FIG. 4. The FLL order parameter in a twinned 11 and in our 
untwinned sample in an applied field of 4 T. In the untwinned 
sample, the melting transition is step-like with a width of <1 
K. In the twinned sample in contrast, the order parameter 
decreases continuously over more than 5 K below the melting 
temperature. Close to T m , the decrease in the twinned sam- 
ple is linear, as would be expected for an order parameter in 
mean field theory. 
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FIG. 5. The /xSR lineshapes above and below the melting 
temperature in an applied field of 0.6 T. At low temperatures, 
the field distribution is asymmetric towards high fields, as is 
expected for a lattice of straight vortex lines due to the high 
local fields in the vortex cores. Above T m however, the line- 
shape is considerably narrowed due to thermal fluctuations. 
The expected symmetric shape is changed to a weighting to- 
wards low fields because of a macroscopic distribution of fields 
over the whole of the sample (see the text). 
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FIG. 6. The temperature dependence of the mean inter- 
nal field (B) and the lineshape asymmetry parameter a at a 
field of 0.3 T. Even at such a low field, there is a sharp melt- 
ing transition as evidenced by the change in sign of a. This 
takes place at a temperature below T c as determined from 
the onset of diamagnetic behaviour in (B). This change in a 
corresponds to that in lineshape displayed in Fig. 5 consisting 
in the truncation of the high field tail due to a smearing of the 
core fields because of the thermal fluctuations of the pancake 
vortices. 
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FIG. 7. The melting line as determined with the field paral- 
lel to the c-direction. Filled circles represent the irreversibility 
line as determined from the onset of irreversible magnetisation 
using a VSM. Open circles are determined by SANS from the 
jump in the FLL order parameter. The squares correspond to 
/iSR measurements, where the change in sign of a was taken 
as an indication of the melting transition. The solid line is a 
fit to all of the data of a melting line of the form B m = Bo 
(1 - T/T c ) n with Bo = 100(10) T and n = 1.30(5). This is in 
agreement with the expectations for a melting line given the 
temperature dependence of A from Fig. 2 (see the text). 
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FIG. 8. The melting line as determined by SANS for differ- 
ent angles with respect to the c-axis. The lines correspond to 
the expectations from the anisotropic scaling theory of Blat- 
ter et al. 9 given the melting line parallel to c (see Fig. 7) and 
an anisotropy of 7 = 4 (see Ref. 18 ). A determination of 7 
from the melting lines is not accurate enough, as the studied 
angles only extend to 50° . 
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FIG. 9. The FLL order parameter in a field of 5 T at 
three different angles to the c-direction. The temperatures 
are scaled to the respective melting temperatures to give a 
unified picture. Within the errors, the transitions are the 
same, which is consistent with the finding of a constant en- 
tropy change as a function of angle 8 . 
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